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ABSTRACT 

In this paper we compute rest-frame extinctions for the aftergfows of a sample of 
Swift 7-ray bursts complete in redshift. The selection criteria of the sample are based 
on observational high-energy parameters of the prompt emission and therefore our 
sample should not be biased against dusty sight-lines. It is therefore expected that our 
inferences hold for the general population of 7-ray bursts. Our main result is that the 
optical/near-infrared extinction of 7-ray burst afterglows in our sample does not follow 
a single distribution. 87% of the events are absorbed by less than 2 mag, and 50% suffer 
from less than 0.3-0.4 mag extinction. The remaining 13% of the afterglows are highly 
absorbed. The true percentage of 7-ray burst afterglows showing high absorption could 
be even higher since a fair fraction of the events without reliable redshift measurement 
are probably part of this class. These events may be due to highly dusty molecular 
clouds/star forming regions associated with the 7-ray burst progenitor or along the 
afterglow line of sight, and/or to massive dusty host galaxies. 

No clear evolution in the dust extinction properties is evident within the redshift 
range of our sample, although the largest extinctions are at z ~ 1.5 — 2, close to 
the expected peak of the star formation rate. Those events classified as dark are 
characterized, on average, by a higher extinction than typical events in the sample. A 
correlation between optical/near-infrarcd extinction and hydrogen-equivalent column 
density based on X-ray studies is shown although the observed Nu appears to be well 
in excess compared to those observed in the Local Group. Dust extinction does not 
seem to correlate with GRB energetics or luminosity. 
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1 INTRODUCTION 

The study of the environments of long-duration 7-ray bursts 
(GRBs) is a subject of growing importance for its implica- 
tions on many different research areas from GRB physics to 
host galaxy chemical evolution. The cosmological nature of 
these events makes them even more interesting, allowing to 
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face the fundamental problem of determining how circum- 
burst environmental parameters vary with the age of the 
universe (see Gehrels et al. 2009, for a review on this field). 

One possible way to obtain information about GRB en- 
vironments along the line-of-sight (LOS) is to study their 
optical/near infrared (NIR) spectral energy distribution 
(SED). This allows researchers to derive information about 
the intrinsic spectrum in this spectral range and to study 
the rest-frame extinction curve, a precious source of infor- 
mation about dust formation (and possibly destruction) in 
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high-redshift environments (e.g. Waxman & Draine 2000; 
Pruchter et al. 2001; Perna & Lazzati 2002; Draine & Hao 
2002; Pcrna ot al. 2003; Stratta et al. 2004; Chen et al. 2006b; 
Lu et al. 2011). 

One of the first attempts to derive statistical informa- 
tion was carried out by Schady et al. (2007) using Swift 
data and their main results are still in agreement with our 
present knowledge: in most cases an extinction curve typical 
of the Small Magellanic Cloud (SMC) environment provides 
a good fit to the data and only in a few cases the extinc- 
tion bump at about 2175 A, typical of the Milky Way (MW) 
and the Large Magellanic Cloud (LMC) environments (e.g. 
Pei 1992) was singled out. Schady et al. (2007) gave also 
support to the idea that afterglows undetected by Swift- 
UVOT\ in several cases classifiable as "dark" bursts (see 
e.g. Melandri et al. 2012, and references therein for a discus- 
sion), were those with high visual extinction (Lazzati et al. 
2002). Schady et al. (2010), with a larger sample of events, 
studied also the effect of soft X-ray absorption, which is usu- 
ally quantified with the hydrogen-equivalent column density 
A'"h. They found that the ratio Nh/Av, with a large scatter, 
is significantly larger than what measured in the Magellanic 
Clouds or the MW. 

These studies however suffered from the limited spectral 
response and sensitivity of Swift-WOT and the high level of 
inliomogerieity of the available observational data for GRB 
afterglows (Kann et al. 2010, 2011). The situation improved 
in recent years thanks to the growing attitude of teams in- 
volved in GRB follow-up to share their data for a better 
coverage of a specific event (e.g. Covino et al. 2008; Jin et 
al. 2013), and in particular for the advent of advanced in- 
struments such as GROND^ equipping the 2.2 m MPI/ESO 
telescope at La Silla (Chile). With its simultaneous 7-band 
imager, from the optical to the NIR, GROND allows to ob- 
tain homogeneous and uniformly calibrated data allowing to 
deal with afterglow SED determination with unprecedented 
reliability (e.g. Greiner et al. 2011). 

The search for highly extinguished GRB afterglows re- 
quires high-quality data of adequate spectral range, and in- 
deed studies devoted to this subject (e.g. Perley et al. 2009; 
Kriihler et al. 2011) could identify several events with mod- 
erate {Ay ~ 1 mag) to high {Ay > 2 mag) rest-frame extinc- 
tion. Greiner et al. (2011), further addressing the problem of 
the nature of dark GRBs, confirmed that for a large fraction 
of them (about 3/4) extinction is responsible for the faint 
optical fluxes, although still about up to 1/4 of the studied 
events are consistent with being high-redshift GRBs, with 
the optical flux depressed by intergalactic neutral hydrogen 
absorption. 

Even better results could in principle be obtained 
through optical/NIR spectroscopy, although the modest 
covered spectral range might sometimes negatively affect 
SED studies. Zafar et al. (2011) studied a sample of 41 
GRBs mainly observed with FORS'* equipping the VLT. Af- 
ter modeling the SEDs from optical to X-rays by means of 
broken power-laws, they remarkably found that in about half 
of the studied events the SEDs require a break between the 



two bands consistent with a synchrotron origin, as required 
by the so-called GRB standard model (Piran 2004). In the 
remaining cases (apart from one outlier) a single power-law 
provides a satisfactory fit. In 63% of the cases a SMC ex- 
tinction curve is preferred, and the 2175 A extinction bump 
is present just in 7% of their sample. About a quarter of 
events is finally consistent with no absorption. 

Any sample so far considered is in any case likely biased 
toward optically bright events. Therefore it is not probably 
totally surprising to see such a high fraction of virtually 
unreddened afterglows. This has clear impact on any in- 
ference about the nature of long GRB progenitors, usually 
supposed to be massive stars still in their star forming re- 
gions (Schulze et al. 2011). A likely improvement of SED 
modeling by means of spectral data will be possible when a 
large sample of afterglow observations, carried out with in- 
struments such as the ESO X-shooter^ (Vernet et al. 2011), 
capable to cover with low-to-medium resolution the whole 
optical/NIR range in one shot, will be available (see, e.g., 
D'Elia et al. 2010; Wiersema 2011). 

From a different point of view, a further step toward 
the characterization of extinction properties of GRB LOSs 
should rely on statistical considerations based on a com- 
plete sample of events according to some selection criterion. 
To this aim we are studying in this paper the SEDs and 
extinction properties of a sample of 58 GRB afterglows (the 
BAT6 sample) described in Salvaterra et al. (2012). All 58 
GRBs in our sample have been selected to have the 1-s peak 
photon fiux P ^ 2.6phs~^ cm~^ (see Salvaterra et al. 2012, 
for further details) as detected by Swift-JiAT. The resulting 
sample has a redshift completeness level of ~ 95% (~ 97% 
of the bursts have a constrained redshift). Several papers 
discussing various features of interest for GRB astrophysics 
for this sample have been delivered so far: X-ray spectra and 
absorbing columns are studied in Campana et al. (2012), the 
percentage of dark burst population is derived in Melandri 
et al. (2012), spectral-energy correlations are discussed in 
Nava et al. (2012) and Ghirlanda et al. (2012), and gen- 
eral prompt/afterglow brightness correlations are derived in 
D'Avanzo et al. (2012). 

The selection criterion based on high-energy brightness 
of the prompt emission avoids the introduction of biases re- 
lated to the optical afterglow detection and, often, for later 
host galaxy identification. Yet, in general, a bright prompt 
emission implies a bright afterglow at any wavelength, al- 
though with a large scatter (D'Avanzo et al. 2012). There- 
fore in a large fraction of cases the optical data available 
are adequate for the analysis. The very high level of red- 
shift completeness makes it possible for the first time to 
obtain representative statistics of the extinction properties 
of GRB LOSs from nearby {z ~ 0.1) to far cosmological 
events {z ~ 5.5). 

In Sect. 2 we discuss the methodologies applied to the 
analysis of the events in our sample. In Sect. 3 we report 
about our main results and in Sect. 4 general conclusions 
are drawn. Details about the analysis for each event are 
discussed in Appendix A. 
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2 DATA AND METHODOLOGY 

In order to compute the SED for as many events as possible 
included in our sample, we collected all the available data 
including those published only in GCN* short communica- 
tions. In a few cases, we had access to still unpublished data 
or derived a new calibration (see Appendix A). 

The SED in the optical/NIR range is modeled as 
a power-law reddened by rest-frame extinction, f„ oc 
i^^'^e^^'^'^^^'" , where is the observed frequency, z is the 
source redshift, t{i') is the computed optical depth and (3 
is the spectral index. We applied three possible extinction 
curves, namely those typical for the MW, LMC and SMC, 
as modeled by Pei (1992). MW extinction in our Galaxy is 
also considered (Sclilcgcl ct al. 1998). Wc did not try a more 
complex modeling of the extinction curve (see, e.g. Calzetti 
et al. 1994; Maiohno et al. 2004; Gallerani et al. 2010; Liang 
& Li 2010) since the inhomogeneity of the available data 
could easily introduce spurious results not justifying the in- 
crease in complexity in our statistical analysis. In addition, 
a reliable analysis of a grey extinction curves would be pos- 
sible only for those events with optical/NIR and soft X-rays 
lying on the same spectral segment, due to the degeneracy 
between optical/NIR spectral normalization and the break 
frequency location between these two spectral regions. This 
would have limited the sample size of our analysis to just 
a small fraction of the available events. Temporal decays 
are modeled as power-laws, f{t) oc (t — to) where t ~ to 
is the time after the burst and a the temporal index. In 
order not to artificially neglect the possible covariance be- 
tween the spectral index and temporal decay when data are 
poorly sampled, temporal and spectral data are simultane- 
ously fit. minimization is performed by using the downhill 
(Nelder-Mead) simplex algorithm as coded in the python^ 
scipy. optimize'^ library, v. 0.10.0. Error search is carried 
out following Cash (1976) with /3 and Eb-v as parameters 
of interest. Throughout this paper the reported uncertain- 
tics arc at la level with the exception of hydrogen-equivalent 
column densities and X-ray spectral slopes, which are shown 
with 90% errors. Upper limits are at 95% confidence level. 
Lower limits have been derived by making use of the usually 
poor information available in the literature, as for instance 
rough upper limits published in GCNs. So, for lower limits, 
it is not possible to give a confidence level. They are anyway 
derived assuming the less demanding, in terms of required 
optical extinction, extrapolation of X-ray spectra to the op- 
tical band, i.e. introducing a spectral break just below the 
Swift-XRT range. 

Our goal is to derive a reliable SED and extinction eval- 
uation, therefore we did not try to model the whole light- 
curve of the analyzed events. Wc instead considered only the 
time interval required to firmly constrain the temporal index 
around the epoch when the most reliable photometric infor- 
mation is available. In several cases, not optimal calibration 
and/or inhomogeneity in the available data required to add 
a systematic calibration error in quadrature. Wc considered 
acceptable a fit if its null probability is better than 10%. 
The spectral slope was derived both considering only opti- 

^ http;//gcn. gsfc.nasa.gov/ 
^ http://www.python.org 
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cal/NIR data and adding a statical prior from the analysis 
of roughly simultaneous X-ray data. This is an alternative 
procedure compared to the simultaneous fit of X-ray and 
optical data, since in this latter case the solution is usu- 
ally dominated by the much more numerous X-ray points, 
making fit results less sensitive to the information provided 
by the optical/NIR data. We imposed that /3o = /3x — 0.5 
or /3o = fix including the 90% error on the X-ray spectral 
slope. This is essentially equivalent to assume that the af- 
terglow emission is due to synchrotron radiation and the 
cooling frequency is either outside or in between the X-ray 
and optical/NIR ranges (see, e.g., Sari et al. 1998; Zhang & 
Meszaros 2004, and references therein). However, this choice 
does not imply that the afterglow evolution can be modeled 
in the framework of the cosmological fireball model (see also 
Uhm & Zhang 2013). Results for all extinction recipes are 
reported in Table AlO (available only in electronic format) 
and details about each event fit and data selection are re- 
ported in Appendix A. In case more solutions are statisti- 
cally acceptable we always chose the one obtained assuming 
a SMC extinction curve, for an easier comparison with most 
of the previous studies of GRB afterglow optical/NIR SEDs. 
Multiple acceptable solutions are typically possible only in 
case of low total extinction, when results based on different 
extinction recipes tend to be comparable. 

For the analysis of the events in our sample we adopted 
the redshifts reported in Salvaterra et al. (2012). X-ray 
spectral indices are taken from Melandri et al. (2012) and 
hydrogen-equivalent column densities from Campana et al. 
(2012). Energetics and luminosities for the prompt emission 
are taken from Nava ct al. (2012). 

The best fit results, which are the base for our further 
analysis, are reported in Table 1, together with redshifts, 
hydrogen-equivalent column densities and spectral slopes. 



3 RESULTS AND DISCUSSION 

In Table 1 wc report results obtained leaving the spectral 
slope of the optical/NIR data completely free or constrain- 
ing this parameter to be in agreement with the results ob- 
tained by the analysis of simultaneous X-ray data. However, 
only in a minority of cases optical data alone allowed us to 
obtain useful estimates of the rcst-framc extinction. In those 
cases results from both analyses are in reasonable agreement 
(Fig. 1). In the rest of the paper we will use only data ob- 
tained using the X-ray information. 

Our sample consists of 58 events, covering the redshift 
range from ~ 0.1 to ~ 5.5. Rest-frame extinctions or lim- 
its on it could be estimated for 53 events (~ 91%). For the 
remaining cases no information could be obtained due to 
the poor quality of optical data, X-ray data, or the lack of 
a redshift. In at least 8 cases 15%) an extinction curve 
presenting a bump at ~ 2175 A, i.e. LMC or MW, is re- 
quired. Zafar et al. (2011), through GRB afterglow spec- 
troscopy, measured this percentage to be about 7%, although 
in their sample high extinction events would likely be under- 
represented for observational biases. Broad-band photome- 
try is clearly not the best tool to single out the extinction 
bump, in particular for low total extinction, and the true 
percentage is likely to be affected by the application of pos- 
sible different extinction recipes (e.g. Calzetti et al. 1994; 
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Figure 1. Comparison between Ay derived witii and without 
assuming an X-ray prior. Only data for which at least one of the 
plotted quantities is not a limit are shown. Here and for later 
figures, blue (circle) symbols are actual measurements and red 
(triangles) symbols are upper or lower limits. 



Maiolino et al. 2004; Liang & Li 2010). In a large fraction 
of cases, a chromatic extinction law is definitely required by 
the data. 

The distribution of rest-frame extinctions is shown in 
Fig. 2 and resembles analogous histograms discussed by 
Kann et al. (2010) and Greiner et al. (2011). However, since 
our sample is selected with purely observational criteria of 
the prompt 7-ray emission, and it is highly complete in 
redshift, is likely to provide an unbiased view of the true 
"parent" distribution of rest-frame extinction for cosmolog- 
ical GRBs. The distribution seems to be bimodal, with a 
smooth distribution clustered at virtually zero extinction 
(including the upper limits) and several events (including 
the lower limits) at higher, roughly Ay ^ 2 mag, extinc- 
tion. 87% of the GRB afterglows (46/53) have an extinction 
smaller than about 2 mag, and 50% smaller than 0.3-0.4 mag. 
The distribution is indeed strongly peaked at low extinction. 
The remaining 13% of events are on the contrary suffering 
from high extinction. They appear to follow, within the lim- 
its of the size of the sample, a different distribution being 
essentially inconsistent with a simple extrapolation of low- 
extinction events at more than 3cr level, if we model the 
low-extinction peak with a Gaussian or a Poissonian func- 
tion. The percentage of events with large extinction could 
be even higher since the few events with no redshift mea- 
surements in our sample are likely part of this class. The 
true percentage of GRB afterglows showing high extinction 
is therefore probably between 15-20%, well in the range of 
the estimated percentage of dark bursts derived by means of 
different criteria and in different samples(de Pasquale et al. 
2003; Jakobsson et al. 2004; Rol et al. 2005; van der Horst 
et al. 2009; Greiner et al. 2011; Melandri et al. 2012). 

The shape of the low-extinction event distribution is 
roughly consistent with what could be expected in typical 
Galactic molecular clouds, although the amount of extinc- 
tion is substantially lower than expected (e.g. Reichart & 
Price 2002; Vergani et al. 2004). The 15-20% of events with 
high extinction could be part of a different population of 



Figure 3. Extinction versus redshift. No clear relation appears to 
be present in the redshift range covered by our sample, although 
the largest extinctions are at z ^ 1.5 — 2, close to the expected 
peak of the star formation rate (z ~ 2 — 3, Hopkins & Beacom 
2006). Bigger symbols (gold: measurements, green: limits) show 
data for dark bursts in our sample according to Melandri et al. 
(2012). 



GRBs born in very dusty molecular clouds (Prochaska et 
al. 2009; Perley et al. 2009) and/or hosted by dusty, more 
massive, galaxies than the average of the host GRB popu- 
lation (e.g. Kriihler et al. 2011; Rossi et al. 2012; Perley et 
al. 2013). 

The derived rest-frame extinctions do not show a clear 
redshift dependence at least up to 2 ~ 4, where we have 
a sufflcient number of events (Fig. 3). The most stringent 
upper limits are at low redshift, but this is probably an 
observational bias since the best optical data are found 
for those events. The largest extinctions are however at 
2; ~ 1.5 — 2, close to the expected peak of the star forma- 
tion rate (z ^ 2 — 3, Hopkins & Beacom 2006). This might 
have interesting consequences although the observed extinc- 
tion in GRB afterglows is related to their LOSs and their 
physical properties do not necessarily correlate with global 
properties of their host galaxies. In the same plot we show 
the position in the Av vs z plane of the dark bursts in our 
sample (Melandri et al. 2012). It is clear that, as expected, 
most dark bursts are characterized by a substantial optical 
extinction, in agreement with previous studies (e.g. Greiner 
et al. 2011; Kriihler et al. 2011). An interesting exception is 
provided by GRB 081222 that is classified as dark (Melandri 
et al. 2012). Yet, it seems to be characterized by a very low 
dust extinction although the classification of this event as 
dark is questionable, as discussed in Melandri et al. (2012). 
On the contrary, GRB 060306 is characterized by a high ex- 
tinction and it was not classified as dark. For this event 
only upper limits in the NIR are available (see Appendix A), 
and therefore even its classification is very uncertain. In any 
case, the classification of an event as dark always requires a 
careful comparison between optical and X-ray emission, and 
therefore the amount of optical extinction is only one of the 
important factors involved in the classification. 

The relation between Av and A^h derived by X-ray anal- 
ysis has been widely discussed in the literature (e.g. Stratta 
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Figure 2. Histogram of the computed rest-frame Ay- Upper and lower limits are also reported. 87% of the GRB afterglows suffer 
from an extinction lower than about 2 mag, and 50% lower than 0.3-0.4 mag. About 13% of the population is characterized by a large 
extinction. 



et al. 2004; Watson et al. 2007; Schady et al. 2010; Zafar 
et al. 2011; Watson et al. 2013). Taking into account the 
metallicity effect on the observed hydrogen-equivalent col- 
umn densities there is a general consensus about a low dust- 
to-gas ratio (~ 10%) compared to typical Local Group (LG) 
values. Our sample (Fig. 4) indeed singles out a well defined 
trend with the highest rest-frame extinctions derived for the 
highest hydrogen-equivalent absorptions, roughly following 
the relation Nu/Av ~ 1.6 x 10^^ cm"^ mag"^ The typical 
LG dust-to gas ratio (Welty et al. 2012) for Solar metallicity 
is also shown. Nonetheless, up to Ay ~ 1 mag the trend is 
not evident and the there is a scattered distribution. The 
hint for two populations found from the Ay distribution 
(Fig. 2) might be found also in the Nn/Ay plan. The most 
absorbed afterglows (with Ay ^ 1 mag) are all characterized 
by the highest A^h values, i.e. larger than lO^'^cm"^. 

Hydrogen-equivalent column densities derived through 
analysis of GRB afterglow X-ray data are actually due to 
photoelectric absorption by inner shells of elements as O, Si, 
S, Fe and He (see Watson et al. 2013, and references therein). 
The relation with hydrogen is obtained typically assuming 
a gas composition in Solar proportion and Solar metallicity. 
The comparison with dust absorption depends also on the 
dust-to-metal ratio (e.g. Savaglio et al. 2003; Zafar & Watson 
2013). None of these parameters can, in general, be robustly 
determined by observations or predicted theoretically for the 
whole sample, making any firm inference about the relation 
between the X-ray A^h and Ay still difficult. Dust sublima- 
tion by the GRB emission (e.g. Waxman & Draine 2000) is 
usually invoked as a possible explanation for the low dust 
content along GRB LOSs, while the effect on the measured 
A^H can be less important since (mildly) ionized metals can 
still produce absorption in low-resolution soft X-ray spectra 
(e.g. Perna et al. 2003). However, as shown in Fig. 5, the 
energetics or luminosity of the GRB is not correlated with 
the amount of dust absorption, as it could be expected in 
case of sizable dust sublimation. 

Watson et al. (2013) reported an evolution in redshift 
for the Nii/Ay relation, with higher values at high redsfhits. 




10.00 



Figure 4. Relation between Ay and A^h derived by X-ray analy- 
sis. Typical dust-to-gas ratio for LG environments is also shown. 
Only data for which at least one of the plotted quantities is not 
a limit have been shown. No correction for metallicity different 
from Solar was applied to the observed A^h- 



Based on our sample (Fig. 6) we can say that this is a direct 
consequence of lack of low A'^h values at high redshift (see 
Campana et al. 2010, 2012), whereas Ay does not evolve 
with redshift (see Fig. 3). The physical reason of the prefer- 
ence for high A'^H values at high redshifts and its significance 
have still to be determined. 

In Campana et al. (2012) this phenomenon is attributed 
to the increasing absorption by intervening systems in higher 
redshift GRBs and therefore it would not be related to 
the progenitor environment (see, however, Watson et al. 
2013 and Starling et al. 2013). The effect is important. 
The A'h contribution due to intervening systems, derived 
by a rough estimate from, e.g. Fig. 2 in Campana et al. 
(2010), at 2 ~ 3 is ~ 6 X 10^° cm"^ and at 2 ~ 6 rises to 
The subtraction of this contribution from 
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Figure 6. Nyi/Ay ratio versus redshift. A lack of low hydrogen- 
equivalent column densities at high redshifts mimics an increase 
of Nn/Ay (see Campana et al. 2010, 2012). 
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Figure 5. Extinction versus E^^a (upper panel) and Liso (lower 
panel). The amount of extinction does not seem to be affected by 
the energetics or luminosity of the GRB. 



the measured hydrogen-equivalent column densities does not 
however make the relation with the optical extinction more 
defined. The intrinsic uncertainties in deriving this relation 
added to the uncertainties associated to the derived A''h 
make it actually difficult to derive firmer conclusions. 

The hydrogen column density can also be measured 
more directly when Lyman-a is visible in GRB afterglow 
spectra. This practically requires that events be at a red- 
shift sufficiently high to have the Lyman-a transition in the 
useful spectral range of ground-based instruments, i.e. from 
the ultraviolet atmospheric cutoff at ~ 3500 A upwards. The 
redshift limit therefore is approximately z > 2. This limits 




xlO 



Figure 7. Relation between A'^h from X-ray spectral analysis and 
Ajji derived by measuring Lyman-«. Blue and red symbols are as 
in Fig. 1. Gold and green symbols are measurements and upper 
limits after subtracting from the X-ray hydrogen-equivalent col- 
umn density a contribution from intervening systems (see Cam- 
pana et al. 2012). 



the number of events for which Lyman-a has been mea- 
sured. In our sample there are just 11 such events. Data 
are from Fynbo et al. (2009, and references therein) but 
for GRB081203A (Kuin et al. 2009). In Fig. 7 we show the 
relation between iVn based on Lyman-a observation (A''hi) 
and A''h based on X-ray analysis. If both A^hi and the gas 
responsible for the X-ray absorption come from the diffuse 
interstellar medium of the host galaxy, the two values should 
be correlated. For some of the events in our sample, rather 
surprisingly, the agreement is satisfactory. Possibly, for the 
sub-sample of events with both measurements, metallicity 
effects on the observed X-ray based A^h are not dominant. 
Nonetheless, as already discussed, the X-ray values are cal- 
culated using Solar metallicity, which might not correspond 
to the host metallicity, and past studies using larger sam- 
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pies found an excess of the X-ray hydrogen-equivalent col- 
umn density compared to A'hi • It is interesting to note that 
if wo again subtract to the X-ray values the mean contri- 
bution in the observed column density resulting from in- 
tervening systems as a function of redshift as empirically 
calculated by Campana et al. (2012), the relation between 
/Vh and A'hi does become tighter (see Fig. 7) , as expected 
if the gas responsible for both absorptions comes only from 
the GRB host. This can be a further hint, within the limits 
of the small available sample, for a possible contribution of 
the intervening systems to the hydrogen-equivalent column 
density determined from the X-ray. 



4 CONCLUSIONS 

In this paper we have computed rest-frame extinctions for 
a sample complete in redshift of GRBs (Salvaterra et al. 
2012). The selection criteria of the sample rely only on ob- 
servational high-energy parameters of the prompt emission 
and therefore our sample is not biased against dusty sight- 
lines. It is expected that our inferences hold for the general 
population of GRBs. 

Our main result is that the extinction suffered by the 
analyzed GRBs does not seem to follow a single distribu- 
tion. 87% of the events are absorbed by less than about 2 
magnitudes, and 50% only suffer from less than 0.3-0.4 mag 
extinction. 13% of the events are instead more absorbed, and 
this population of GRBs may be due to highly dusty molec- 
ular clouds/star forming regions associated with the 7-ray 
burst progenitor or along the afterglow line of sight, and/or 
to massive dusty host galaxies. The true percentage of GRB 
afterglows showing high absorption is probably higher, since 
most of the events without reliable redshift measurements 
are likely part of this class. 

No clear evolution in the dust extinction properties is 
evident within the redshift range of our sample (~ 0.1 — 5.5), 
although the largest extinctions are at a ~ 1.5 — 2, close to 
the expected peak of the star formation rate (2 ~ 2 — 3, 
Hopkins & Beacom 2006). Those GRBs classified as dark 
(Melandri et al. 2012) show higher than average extinction, 
confirming past claims about the nature of under luminous 
optical afterglows. There is a well defined relation between 
optical/NIR extinction and X-ray hydrogen-equivalent col- 
umn density, although with a gas-to-dust ratio well above 
that observed in LG environments. Dust extinction does not 
seem to correlate with GRB energetics or luminosity, as it 
could be expected if dust properties are affected by the pho- 
ton flux from the high-energy event. 
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Table 1. Spectral and absorption parameters for GRBs in our sample. For each event we list the rcdshift, the X-ray spectral index, the 
X-ray hydrogen-equivalent absorbing column, the optical/NIR spectral index and the amount of rest-frame reddening in the visual band. 
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This figure is reported in Greiner et al. (2011), although the related photometry has not been published. 
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APPENDIX A: GRBS ANALYSED IN OUR 
SAMPLE 

In this section we will frequently refer to various photometric 
bands. With capital letters we refer to the Johnson-Morgan- 
Cousin bands: U,B,V,Rc,Ic- With primed small letters 



we refer to the SDSS photometric system u',g',r',i',z'^. 
Swift-WOT bands are again referred with lower-case let- 
ters uvw2, uv'm2, uvwl, u, b, v, wh. In a few cases systematic 
errors to bring the fit quality to a level adequate for error 
modeling were introduced. Which data were selected from 
the literature and specific comments are reported below for 
each analyzed event. 

GRB 050318 Swift-VYOT u, b and v (AB) magnitudes, 
reported in Still et al. (2005), and an Rc band measurement, 

reported by Bcrgcr ct al. (2005), arc available. Calibrations 
of these two datasets do not appear to be in agreement, and 
therefore we neglected the Rc point. 

GRB 050401 Wc collected data from Rykoff et al. (2005); 
de Pasquale et al. (2006); Watson et al. (2006); Kamble et 
al. (2009) in the V,Rc,IJc,J,H and K bands. However, 
according to the authors, R band data in de Pasquale et al. 

(2006) arc poorly calibrated and have not been used in the 
analyses. In addition R band data in Watson et al. (2006) 
and Kamble et al. (2009) are in strong disagreement. Given 
the more detailed calibration process reported by Kamble et 
al. (2009) we chose to use these data. 

GRB050416A We collected data from Holland et al. 

(2007) ; Soderberg et al. (2007) and Price et al. (2005) in 
the Swift-XJYOT uvm2,u,b,v and Rc,Ic,z and K bands. 
From this dataset we did not consider the GCN data due to 
the very preliminary calibration and one / point in Holland 
et al. (2007) at about 1150 min that was clearly discrepant. 
We also removed SwiftrUVOT data with error larger than 
0.4 mag. 

GRB 050525A This event was observed intensively. We 

collected Swift-VYOT data in the uvw2,wl.u,b,v bands 
from Blustin et al. (2006). We neglected all data with cali- 
bration based only on preliminary photometry. 

GRB 050802 We collected data from Gates et al. (2007a) 
in the Swift-VYOT uvwl, u, b, v bands and from Table Al in 
the u' , B, V, Rc and 7c bands. In Testa et al. (2005) prelimi- 
nary TNG^ photometry was reported. Here we re-calibrated 
all the available data following the field calibration by Hen- 
den et al. (2005) for EVRch and SDSS** photometry for 
u' . We neglected SwiftAJYOT data obtained after about 
8 hr after the GRB and magnitudes with errors larger than 
0.4 mag. In any case Sm/t-UVOT still shows a considerable 
scatter, likely intrinsic, and to obtain a satisfactory fit we 
had to introduce a 15% systematic error to be added in 
quadrature. 

GRB 050922C We collected data from Gates et al. (2009) 
in the SwiftrUYOT u,b,v bands. In Table A2 we report 

VLT'"* photometry in B,V and Rc bands calibrated by 
Landolt standard field stars and preliminarily published in 



^ http://www.sdss.org/ 

http://www.tng.iac.es/ 
* http:/ /skyserver.sdss3.org/dr8 /en/ 
^ http://www.eso.org/ 
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Table Al. TNG data for GRB 050802. u' magnitude is in the 

AB system, BVRcIc arc in tlic Vega system. 



t-to 


Exposure 


Magnitude 


Band 


Telescope 


(min) 


(s) 








740.0 


720 


22.73 ±0.15 


a' 


TNG 


727.5 


240 


21.81 ±0.04 


B 


TNG 


722.0 


120 


21.37 ±0.05 


V 


TNG 


2135.0 


600 


22.85 ±0.07 


V 


TNG 


716.5 


120 


21.07 ±0.05 


Rc 


TNG 


709.5 


240 


20.44 ±0.06 


Ic 


TNG 


2197.0 


360 


21.90 ±0.15 


Ic 


TNG 



Table A2. TNG data for GRB050922C. Magnitudes are in tlie 
Vega system. 



t-to 

(min) 


Exposure 

(s) 


Magnitude 


Band 


Telescope 


1700.0 


180 


22.41 ±0.03 


B 


VLT 


1686.7 


200 


22.00 ±0.02 


V 


VLT 


1692.4 


120 


21.52 ±0.04 


Rc 


VLT 



Covino et al. (2005). In the fit we used only data earlier than 
about 20 min. 

GRB 060206 This event was studied by several groups 
(Wozniak et al. 2006; Monfardini et al. 2006; Stanek et al. 
2007; Curran et al. 2007a). We collected optical data from 
Haislip et al. (2006); Monfardini et al. (2006); Curran et al. 
(2007a); Stanck ct al. (2007) and Gates et al. (2009) in the 
Swift-WOT V and in the r',Rc,Ic bands, and NIR data 
from Alatalo et al. (2006) and Terada et al. (2006) in the 
J,H,K bands. In Table A3 we report Asiago^" and TNG 
photometry in the r' and i' bands calibrated by SDSS^^ 
photometry of field stars, which was preliminarily published 
in Malesani et al. (2006a). The available data come from sev- 
eral different telescopes and in spite of the efforts devoted to 
cross-calibration by the various authors we had to introduce 
a 5% systematic error to be added in quadrature. In addi- 
tion, B and V bands are bluer than Lyman-a at the redshift 
of this event. As shown, e.g., by Aoki et al. (2009) this spec- 
tral range is affected by a dense Lyman-a forest making it 
difficult to model its absorption. We therefore limited our 
analysis to data redder than Lyman-a. 



http://www.pd.astro.it/asiago/ 
http://skyserver.sdss3.org/dr8/en/ 

Table A3. TNG data for GRB 060206. Magnitudes are in the 
AB system. 



t-to 


Exposure 


Magnitude 


Band 


Telescope 


(min) 


(s) 








1158.5 


2000 


19.95 ± 0.06 


r' 


Asiago 


2695.7 


3600 


21.02 ±0.06 


r' 


Asiago 


11534.1 


3000 


23.52 ±0.35 


r' 


TNG 


1229.2 


2700 


19.48 ±0.05 


i' 


Asiago 



Table A4. VLT data for GRB 060814. The host galaxy is likely to 
dominate these measurements. Magnitudes are in the AB system. 



t-to 
(min) 


Exposure 
(s) 


Magnitude 


Band 


Telescope 


50.1 


900 


24.39 ±0.11 


r' 


VLT 


1469.9 


960 


24.71 ± 0.20 


r' 


VLT 


89.3 


360 


24.11 ± 0.22 


i' 


VLT 


1464.4 


1440 


24.56 ± 0.31 


i' 


VLT 



GRB 060210 Data for this event were published by Stanek 
et al. (2007), Rc, Curran et al. (2007b), Rc and Ic, and Cenko 
et al. (2009), Rc,i' and z' bands. Unfortunately, as already 
pointed out by Kann et al. (2010), there is a calibration offset 
between 7c data from Curran et al. (2007b) and i data from 
Cenko et al. (2009). We had therefore to add in quadrature 
a 5% systematic error. 

GRB 060306 No afterglow detection was reported for this 
event. Upper limits are reported by de Pasquale et al. 
(2006b) in the Swift-WOT bands and by Nysewander et 
al. (2006); Chen et al. (2006a) and Lamb et al. (2006) in 
the NIR. Gnly a lower limit for Av could be derived. The 
redshift for this event was questioned by Jakobsson et al. 
(2012); Kriihler et al. (2012), and a solution at lower red- 
shift is reported by Perley et al. (2013). 

GRB 060614 This low-redshift event was widely followed. 
We collected Swift-UYOT uvw2, uvm2, uvwl, u, b, v and wh 
data from Mangano et al. (2007) and B, V, Rc, Ic, J, K data 
from Delia Valle et al. (2006). These two main datasets show 
small calibration inconsistencies and we had to introduce a 
5% systematic error to be added in quadrature. 

GRB 060814 The afterglow was detected in the NIR 
(Levan et al. 2006a; Cenko et al. 2006a; Levan et al. 2006b), 
with its emission already contaminated by the relatively 
bright host galaxy. The same is true for our r' and i' VLT 
observations calibrated by SDSS photometry of field stars, 
which were preliminarily discussed in Malesani et al. (2006b) 
and reported in Table A4. We derived a lower limit on Ay 
assuming that the afterglow cannot be brighter than these 
detection, likely dominated by the host galaxy contribution. 
With reference to Figure 9 in Jakobsson et al. (2012) we have 
tried to isolate the photometry of object A only. 

GRB 060904A No adequately calibrated observations of 
this afterglow have been published, and the redshift is also 
not known. 

GRB 060908 Data for this event are collected by Cenko 
et al. (2009) {g',Rc,i' and z' bands) and Covino ct al. 
(2010) (Sw/t-UVOT u, V and wh and B, V, Rc, Ic, J, H and 
K bands). The total dataset is remarkable although due to 
the many different telescope/instrument combinations we 
had to add in quadrature an 8% systematic error, as already 
discussed in Covino et al. (2010). 

GRB060912A Data for this event are collected 
from Gates et al. (2009) in the Swift-VVOT 
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uvw2, uvm2, uvwl, u, h, v and wh bands and from Deng et 
al. (2009) in the -Rc band. Lyman-a absorption affecting the 
bluest Swift-WOT filter was modeled superposing to the 
power-law SED a Voigt profile with parameters reported in 
Deng et al. (2009). 

GRB 060927 Calibrated data for this event were reported 
by Ruiz-Velasco et al. (2007) in the i',Rc,Ic and K bands. 
The Rc band flux was corrected for Lyman-a absorption 
with parameters obtained by Ruiz-Velasco et al. (2007). 

GRB 061007 Data are collected from MundeU et al. 
(2007) {Rc band) and Schady et al. (2007) {Swift-WOT 
uvwl,u,b,v and wh bands). Cross-calibrations of the 
datasets required to add in quadrature a 3% systematic er- 
ror. 

GRB 061021 Data are collected from Gates et al. (2009) 
in the Swift-UYOT uvw2, uvm2, uvwl, u, b, v and wh bands. 

GRB 061121 Data are collected from Gates et al. (2009) 
in the Swift-UVOT uvwl, u, b, v and wh bands and from 
Kann et al. (2010) in the V, Rc and Ic bands. We removed 
the uvwl band since it is likely heavily affected by Lyman-a 
absorption. 

GRB 061222A Only K band afterglow detection was re- 
ported by Cenko et al. (2006b, c). The available NIR data 
can not constrain the fit parameters and we derived the ex- 
tinction by extrapolating the X-ray spectrum. 

GRB 070306 Data in the H and K bands are collected 
from Jaunsen et al. (2008). 

GRB 070328 No afterglow detection was reported for this 
event and the redshift is not known. Upper limits in Swift- 
UVOT bands are reported by Marshall et al. (2007). 

GRB 070521 No solid afterglow detection was reported for 

this event. The available photometry refers to the relatively 
bright host galaxy. Upper limits from Minezaki et al. (2007) 
have been used to derive a lower limit on Ay- 

GRB 071020 Data are collected from Cenko et al. (2009) 
in the Rc,i' and z' bands and from Bloom et al. (2007) in 
the J and H bands. The light curve of this event shows a 
considerable scatter, and we had to add in quadrature an 
11% systematic error. 

GRB 071112C Data in the Swift-UYOT uvwl, u, b, v and 

wh bands were collected from Gates et al. (2007b). NIR J 
and K bands data wore collected from Uchara et al. (2010) 
and V, Rc and 7c bands were collected from Huang et al. 
(2012). g',r' and 'i' bands observations obtained at Asiago 
and TNG and calibrated following the SDSS are reported in 
Table A5. 

GRB 071117 The afterglow discovery was reported by 
Bloom et al. (2007b, c). The redshift was derived by Jakob- 
sson et al. (2007) when the host galaxy was likely already 
dominating the afterglow light. Data in Table A6, calibrated 



Table A5. Asiago and TNG data for GRB071112C. Magnitudes 

arc in the AB system. 



t-to 


Exposure 


Magnitude 


Band 


Telescope 


(min) 


(s) 








97.7 


600 


20.88 ±0.15 


g' 


Asiago 


119.7 


300 


20.90 ±0.17 


r' 


Asiago 


126.7 


300 


21.01 ±0.28 


i' 


Asiago 


3274.4 


300 


24.04 ± 0.09 


r' 


TNG 


32036.9 


600 


22.36 ±0.12 


i' 


TNG 


125409.9 


600 


23.93 ± 0.20 


i' 


TNG 



Table A6. VLT and Gcmini-S observations of GRB 071117. Mag- 
nitudes are in the Vega system. 



t-to 


Exposure 


Magnitude 


Band 


Telescope 


(min) 


(s) 








558 


300 


23.04 ±0.10 


Rc 


VLT 


574 


900 


23.02 ±0.04 


Rc 


Gemini-S 


694 


540 


23.02 ±0.05 


Rc 


VLT 


2009 


1440 


23.53 ±0.04 




Gemini-S 



by standard stars observed with the VLT, supersede those 
reported in Bloom et al. (2007b,c) and Jakobsson et al. 
(2007). These data and upper limits from Berger & Gonzalez 
(2007) do not allow to constrain the rest frame extinction. 

GRB 080319B This is the famous "naked eye" burst. Data 
in the V were collected from Pandey et al. (2009). Prom 
Bloom et al. (2009) we used data in r' band and in the 
B,V,Rc and Ic bands. Finally, data in the Swift-UYOT 
uvw2,uvm2,uvwl,u,b,v and wh bands and in the Rc,Ic 
bands were obtained from Racusin et al. (2008). We how- 
ever did not consider in the fit the Swift-UYOT uvw2, uvm2 
and uvwl filters for possible Lyman-a contamination. In ad- 
dition in many cases the reported photometric errors are 
likely neglecting the absolute calibration error judging from 
the unrealistic small quoted uncertainties. Therefore we had 
to add in quadrature a 3% systematic error. 

GRB 080319C Data are coUected from Cenko et al. (2009) 
in the Rc, i' and z' bands. 

GRB080413B Data in the g',r',i' and z' bands, in the 
J, H and Kg bands and in the Swift-UYOT v band were 
collected from Filgas et al. (2011a). We added in quadra- 
ture a 3% systematic error to compensate for the very low 
photometric errors of the available data possibly neglecting 
absolute calibration uncertainties. 

GRB 080430 The only reliably calibrated data pub- 
lished so far for this event are in the Swift-UYOT 
uvw2, uvm2, uvwl, u, b, v and wh bands from Landsman et 
al. (2008). 

GRB 080602 No afterglow detection has been published 
for this event. In addition the available XRT data do not 
allow a meaningful comparison with available optical upper 
limits as in Malesani et al. (2008a). 
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Table A7. REM data for GRB 080916A. Magnitudes are in the 

Vega system. 



t — to Exposure Magnitude Band Telescope 
(min) (s) 



1.27 


50 


13.78 ±0.06 


H 


REM 


2.59 


50 


14.49 ±0.11 


H 


REM 


3.92 


50 


14.12 ±0.08 


H 


REM 


5.26 


50 


14.13 ±0.08 


H 


REM 


6.41 


25 


14.37 ±0.14 


H 


REM 


16.25 


100 


15.14 ±0.15 


H 


REM 



GRB 080603B Data in the Swift-HWOT uvwl, u, h and v 
bands were collected from Kuin et al. (2008). Data in the 
g\r' and i' bands were collected from Jelinek et al. (2012), 
and in the J, H and K bands from Miller et al. (2008) . We 
did not use data bluer than the Lyrnan-a line at this GRB 
redshift due to strong flux suppression (see, e.g. Jelinek et 
al. 2012). 

GRB 080605 Data in the <?', r', i' and z' bands and in the 
J,H and Ka bands were collected from Zafar et al. (2012). 

GRB 080607 Data in the Swift-WOT v band were col- 
lected from Schady ct al. (2008). Data in the V, Rc, Ic, J, J 
and Ks bands, and in the i' and z' bands were collected from 
Perley et al. (2011). We added in quadrature a 5% system- 
atic error to compensate for cross-calibration problems and 
intrinsic short-term variability in the data. 

GRB 080613B No optical afterglow detection has been re- 
ported. 

GRB 080721 Data in the V and Rc, and in the r' and i' 
bands were collected from Starling et al. (2009). We did not 
consider data bluer than Lyman-a. 

GRB 080804 Data in the Swift-VVOT u,b,v and wh 
bands were collected from Kuin et al. (2008b) and data from 
GCN circulars to define the temporal decay only. 

GRB080916A Data in the Swift-WOT u,b,wh bands 
were collected from Gates ct al. (2008); Ziaccpour ot al. 
(2008). H band observations obtained with REM and cal- 
ibrated following 2MASS stars in the field are reported in 
Table A7. 

GRB 081007 Data in the Swift-UYOT u, b, v, wh bands, 
in the B,V,Rc,Ic,H and and in the r' and i' bands were 
collected from Jin et al. (2013). Due to the large number of 

detectors/telescopes involved wc had to add in quadrature 
a 7% systematic error to compensate for cross-calibration 
problems. 

GRB 081121 Data in the Swift-VYOT u,b,v and wh 
bands were collected from Gates et al. (2008b). Data in 
the g' ,r' ,i' ,z' , J,H and Kg were collected from Loew et al. 
(2008). 

GRB081203A Data in the Swift-VYOT 
uvw2,uvwl,u,b,v and wh bands were collected from 



Table A8. REM data for GRB 081222. J, H and Ks magnitudes 

are in the Vega system, r' magnitudes in the AB system. 



t — to Exposure Magnitude Band Telescope 
(min) (s) 



1.37 


50 


11.24 ± 0.09 


H 


REM 


2.69 


50 


12.24 ± 0.10 


H 


REM 


4.02 


50 


12.95 ± 0.07 


H 


REM 


10.01 


30 


15.94 ±0.11 


r' 


REM 


10.67 


30 


16.00 ±0.11 


r' 


REM 


11.32 


30 


16.10 ±0.12 


r' 


REM 


11.98 


30 


16.28 ±0.13 


r' 


REM 


12.61 


150 


14.57 ± 0.07 


U 


REM 


12.63 


30 


16.28 ±0.14 


r' 


REM 


12.29 


30 


16.30 ±0.14 


r' 


REM 


13.94 


30 


16.36 ±0.15 


r' 


REM 


14.59 


30 


16.20 ±0.16 


r' 


REM 


17.27 


30 


16.51 ±0.17 


r' 


REM 


17.92 


30 


16.68 ±0.19 


r' 


REM 


18.58 


30 


16.62 ± 0.18 


r' 


REM 


19.23 


30 


16.73 ± 0.20 


r' 


REM 


19.89 


30 


16.63 ±0.18 


r' 


REM 


28.08 


150 


15.93 ±0.15 


J 


REM 



de Pasquale et al. (2008) and Kuin et al. (2009). Rc data 
were collected from FatkhuUin et al. (2008). We did not 
consider data bluer than Lyman-a. 

GRB 081221 No afterglow detection has been reported for 
this event. The upper limit reported by Malesani et al. 
(2008b) was used to derive a lower limit for Ay . 

GRB 081222 Data in the Swift-VYOT u,b,v and wh 
bands were collected from Brcevold et al. (2008). Data in 
the J,H and bands from Updike ct al. (2008). REM 
photometry is reported in Table A8. These data update and 
substitute those reported in Covino et al. (2008b). NIR data 
were calibrated with 2MASS stars in the field and optical 
data with stars reported in the APASS cataloguc^^. Wc did 
not consider data bluer than Lyman-Q. We had to add in 
quadrature a 7.5% systematic error to compensate the unre- 
alistic low errors reported in some of the available photom- 
etry possibly neglecting absolute calibration uncertainties. 

GRB 090102 Data in the Swift-VYOT u,b,v and wh 
bands were collected from Gurran et al. (2009). Data in the 
g',Rc,i',z',J,H and Ks were collected from Gendre et al. 
(2010). 

GRB 090201 The afterglow detection, likely already sub- 
stantially contaminated by the host galaxy, was reported by 
D'Avanzo ct al. (2009). Data reported in D'Avanzo et al. 
(2009) are superseded by those reported in Table A9. Op- 
tical data were calibrated by observation of standard star 
frames and NIR data following the 2MASS catalogue. The 
redshift for this event is reported in Kriihler et al. (2013b). 

GRB 090424 Data in the Swift-VYOT u,b,v,wh bands 
and in the Rc and Ic were collected from Jin et al. (2013). 

http://www.aavso.org/apass 
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Table A9. VLT data for GRB 090201. Magnitudes are in the 

Vega system. 



t-to 
(min) 


Exposure 
(s) 


Magnitude 


Band 


Telescope 


434.0 


720 


24.12 ±0.08 


Rc 


VLT 


916574 


5000 


24.40 ± 0.09 


Rc 


NTT 


451.0 


480 


23.75 ± 0.11 


Ic 


VLT 


554.0 


1200 


22.05 ± 0.13 


J 


VLT 


580.4 


1200 


21.31 ±0.15 


H 


VLT 


608.0 


1200 


20.19 ±0.14 


Ks 


VLT 



Data in the g' ,r' ,i' , z' , J, H and Ks were collected from 
Olivares et al. (2009). Data in the V, Jc, J and Ks bands 
from Cobb (2009a). Due to the large number of detec- 
tors/telescopes involved we had to add in quadrature a 7.5% 
systematic error to compensate for cross-calibration prob- 
lems. 

GRB 090709 A Data in the z',J,H and Ks bands were 

collected from Ccnko et al. (2010). The rcdshift for this event 
is reported by Perley et al. (2013). Data show a considerable 
scatter, and we had to add in quadrature a 9% systematic 
error. 

GRB 090715B Data in the and Ic bands were collected 
from Haislip et al. (2009), data in the r', i' and z' from Virgili 
et al. (2013, in preparation). 

GRB 090812 Data in the Swift-VVOT u and b bands were 
collected from Schady et al. (2009). Data in the g' ,r' ,i',z' 
and j bands were collected from Updike et al. (2009) and 
Virgili et al. (2013, in preparation). We did not use data 
bluer than Lyman-a. 

GRB 090926B No reliably calibrated data have been pub- 
lished for this event although in Greiner et al. (2011) a SED 
analysis was reported. Upper limits from Xu et al. (2009) 
did not allow to constrain rest frame extinction. 

GRB 091018 Data in the Swift-\]YO'T uvwl,u,b,v,wh 
bands, in the g' ,r' ,i' , z' bands, and in the J, H and K bands 
were collected from Wicrsema et al. (2012). Due to the large 
number of detectors/telescopes involved we had to add in 
quadrature a 2.5% systematic error to compensate for cross- 
calibration problems. 

GRB 091020 Data in the Swift-XJVOT uvwl,u,b,v and 
wh bands were collected from Oates et al. (2009b) and 
Racusin et al. (2009) and in the B, V, Ro and Ic from Kann 
et al. (2009a,b). We did not use data bluer than Lyman-a. 

GRB 091127 Data in the g' , r', i' , z' , J and H bands were 
collected from Filgas et al. (2011b), data in the Ic band from 
Vergani et al. (2011), data in the Swift-VYOT uvwl and u 
bands, and in the B, V, Rc and 7c bands from Cobb (2010). 
Due to the large number of detectors/telescopes involved we 
had to add in quadrature a 2% systematic error to compen- 
sate for cross-calibration problems. 



GRB 091208B Data in the Swift-UYOT u and b bands 
were collected from de Pasquale et al. (2009), data in the 
g' ,r' ,i' , z' , J and H bands were collected from Updike et al. 
(2009b). We however did not use data bluer than Lyman- 
a and collected not-properly calibrated data from GCNs to 
constrain the temporal decay only. 

GRB 100615A No detection is reported for this afterglow. 
Upper limits are discussed in D'Elia & Stratta (2011). Red- 
shift was recently measured by Kriihler et al. (2013a). 

GRB 100621A Data in the g',r', i',z',J,II and Ks bands 
were collected from Kriihler et al. (2011) and data in the 
J,H and Kg were collected from Naito et al. (2010). Data 
from Kriihler et al. (2011) are late enough to be considerably 
affected by the host galaxy brightness, in particular in the 
bluest available bands. We deal with this uncertainty by 
adding in quadrature a 7.5% systematic error. 

GRB100728B Data in the Swift-XJYOT u,b and wh 
bands were collected from Oates et al. (2010), data in the 
Rc band from Volnova et al. (2010) and data in the J, H and 
Ks bands from Olivares et al. (2010). 

GRB110205A Data in the Swift-VYOT u,b,v and wh 
bands, in the g',r',i' and z' bands and in the J,II and 
Ks bands were collected from Cucchiara et al. (2011). Data 
in the U, B, V, Rc, Ic,g' , r' , i' and z' bands from Zheng et al. 
(2012). Data in the B, V, Rc and Ic bands from Gcndrc ot al. 
(2012). Data in the g' , Rc and Ic bands were collected from 
Kuroda et al. (2011a,b). Finally, data in the r' band were 
collected from Urata et al. (2011). To compensate for cross- 
calibration problems in different datasets we added a 2% 
systematic error. We did not use data bluer than Lyman-a. 

GRB 110503 A Data in the Swift-VYOT 

uvwl,uvrn2,uvw2,u,b,v and wh bands were collected 
from Oates et al. (2011) and data in the J and H bands 
were collected from Morgan & Bloom (2011). To compen- 
sate for cross-calibration problems in preliminary datasets 
we added a 7.5% systematic error. We did not use data 
bluer than Lyman-a. 

This paper has been typeset from a T^X/ M^rjX file prepared 
by the author. 
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Table AlO. Results of the analysis performed in this work for the GRBs in the sample. We report the spectral index and the amount 
of rest-frame reddening for each applied extinction recipe; MW, LMC and SMC. Spectral slopes and extinctions are computed leaving 
the spectral slope free or imposing a prior from X-ray data as described in Sect. 2. 



Event 


(min) 


Ext. curve 


/So 


No X-ray prior 

Eb-v xVdof 


/9o 


With X-ray prior 

Eb-v X^/dof 


GRB 050318 


43, 860 


MW 
LMC 
SMC 


4 78+1-" 

9 90+1.14 


-0.19 
34+0-29 

r:o+0-40 
'^•°'^-0.35 


1.00 (4.00/4) 
1.00 (4.00/4) 
1.00 (4.00/4) 


1.02 

1 02+°°° 
89+° " 


0.00 

57+°-l'' 
^•O'-0.15 
9P.+0.06 


16.25 (65.01/4) 
1.47 (5.86/4) 
1.25 (5.02/4) 


GRB 050401 


29, 1400 


MW 
LMC 
SMC 


0.70 
-0.41 

-0 42+°-85 


0.00 

0.20 
16+°-i° 


2.16 (12.96/6) 
1.84 (11.04/6) 
1.43 (8.59/6) 


+0.48 
+0.19 

+o.i9l°;'oo 


0.03 
0.11 

0-07l°;°I 


2.41 (14.43/6) 
2.12 (12.71/6) 
1.64 (9.85/6) 


GRB 050416A 


43, 10000 


MW 
LMC 
SMC 


0.87lgf, 
1 04+0-25 


34+0-35 
'■'■■^^-0.27 

0.14l°o:l^o 
O-OStVol 


0.84 (20.17/24) 
0.78 (18.73/24) 
0.78 (18.60/24) 


o.48t°:?t 

,-, 79+0.00 
^•'^-0.25 
f, 79+0.00 
"■'''-0.17 


o.3ot°:?^ 

r, 1O+0.09 

^•l°-0.03 

l6+°-°5 


0.84 (20.20/24) 
0.79 (18.98/24) 
0.86 (20.73/24) 


GRB 050525A 


2.9, 7.8 


MW 
LMC 
SMC 


1.48 
0.73 

-0.52 


0.00 
0.10 
i2+°-°5 

L'--L^_0.04 


2.16 (38.92/18) 
1.81 (32.63/18) 
1.42 (25.63/18) 


0.73 
0.73 


0.17 
0.10 

10+°°° 


3.02 (54.36/18) 
1.82 (32.84/18) 
1.45 (26.06/18) 


GRB 050802 


86, 2900 


MW 

LMC 
SMC 


70+1 08 


< 0.16(0.01) 

< 0.55(0.13) 

< 0.43(0.13) 


1.71 (10.47/6) 
1.51 (9.05/6) 
1.31 (7.85/6) 


0.43 

43+0-00 


0.22 

'J- 1' -0.03 
Q 10 + 0.03 


4.70 (28.23/6) 
1.53 (9.21/6) 
1.31 (7.85/6) 


GRB 050922C 


3.9, 16 


MW 
LMC 
SMC 


1 91 +0.25 
1-^1_0.65 

0.79l«:^« 


< 0.09(0.01) 

< 0.22(0.04) 

< 0.56(0.10) 


1.20 (19.22/16) 
1.20 (19.21/16) 
1.20 (19.21/16) 


0.811;-° 
o.78i;;- 


'J-'^^-o.oi 
04+"'"^ 
03+»'0i 


1.25 (19.92/16) 
1.20 (19.21/16) 
1.20 (19.21/16) 


GRB 060206 


72, 330 


MW 
LMC 
SMC 


n or;+0.07 
'^■°°-0.33 

02+°-®* 
^■"■^-0.60 


< 0.13(0.00) 
11+0-08 

< 0.09(0.04) 


1.76 (8.79/5) 
1.14 (5.72/5) 
1.39 (6.93/5) 


n oc:+0.07 

f, 97+0.50 
" ' -0.00 
46+° ''l 


< 0.13(0.00) 
0.08l°;°i 

< 0.06(0.04) 


1.76 (8.79/5) 
1.20 (6.01/5) 
1.39 (6.93/5) 


GRB 060210 


8.6, 43 


MW 
LMC 
SMC 


4.47+9-^5 

' —4.86 

4.47i°;I? 


< 1.00(0.00) 

< 0.96(0.00) 


1.06 (22.29. /21) 
1.06 (22.29/21) 
1 Ofi ^"99 9Q/91 


-1 -lo+O.OO 
-^•-^ ■'-(). 10 
1.13+"-™ 

1 iq + 0.00 

-L-J-0_Q j^O 


60+°-^° 
0.281°;°^ 
20+° °^ 


1.24 (26.14/21) 
1.13 (23.77/21) 
1 1 n ^93 01 /9i ^ 


GRB 060306 


96 


MW 
LMC 
SMC 








0.86 

U.OU 

0.86 


> 2.8 

> 2.7 

> 3.0 




GRB 060614 


2300,5300 


MW 
LMC 
SMC 


0.07 
0.03+°-^^ 

— U.ZD 

0.l8ig;f3 


0.33 
f, o,+a.07 

U-51_o.07 
26+° '"' 


2.04 (18.40/9) 
1.39 (12.48/9) 
1.36 (12.28/9) 


0.35 
0.35 

„ or;+0.09 

U.OO_Q QQ 


0.24 

0.22 
P, 99+0.02 
"•^^-0.02 


2.38 (21.43/9) 
1.84 (16.90/9) 
1.52 (13.66/9) 


GRB 060814 


50 


MW 
LMC 
SMC 








0.56 
0.56 
0.56 


> 0.4 

> 0.4 

> 0.4 




GRB 060904A 




MW 
LMC 
SMC 














GRB 060908 


8.6, 290 


MW 
LMC 
SMC 


0.38iro? 
0.38+g;« 

22+°'^^ 

U-^^_0.21 


< 0.01(0.00) 

< 0.04(0.00) 

< 0.07(0.02) 


0.99 (87.80/89) 

0.99 (87.80/89) 
0.97 (86.50/89) 


';fi+° °2 

"■^"^-0.00 

56+° °2 


< 0.01(0.00) 

< 0.01(0.00) 

< 0.01(0.00) 


1.14 (101.60/89) 
1.14 (101.60/89) 
1.14 (101.60/89) 


GRB 060912A 


8.6, 290 


MW 

LMC 
SMC 


-0.64 
25+"-'''^ 

-o.osil:^^ 


< 0.41(0.00) 

< 0.53(0.13) 

< 0.41(0.14) 


0.87 (16.44/19) 
0.83 (15.77/19) 
0.77 (14.71/19) 


90+° °° 

46+°-^* 
^•™-o.oo 

46+°-''* 


< 0.22(0.03) 

< 0.16(0.10) 
07+°'°2 

'-'■'^'-0.06 


0.89 (16.90/19) 
0.83 (15.82/19) 
0.79 (15.07/19) 


GRB 060927 


14, 4500 


MW 
LMC 
SMC 


n 75+0-41 
0.75i?;« 
0.751^3^ 


< 0.17(0.00) 

< 0.41(0.00) 

< 0.47(0.00) 


1.75 (6.98/4) 
1.75 (6.98/4) 
1.75 (6.98/4) 


75+0--11 
^•"'-0.04 
75+0-41 

75+°-4l 


< 0.10(0.00) 

< 0.08(0.00) 

< 0.06(0.00) 


1.75 (6.98/4) 
1.75 (6.98/4) 
1.75 (6.98/4) 
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Table AlO - continued 



GRB 061007 


7.2, 22 


MW 


2.36 


0.13 


1.80 (117.23/65) 


1.10 


0.37 


11.05 (718.29/1 






LMC 






1.19 (77.40/65) 


97+0 13 

-0.02 


oQ+o.oi 

^J-^'5_0.01 


1.20 (77.76/6J 






SMC 


-0.30l°;« 


24+0 03 


1.18 (76.94/65) 


0.95 


0.14 


1.44 (93.84/61 



GRB 061021 


58, 360 


MW 


40+0 09 


< 0.14(0.00) 


0.45 (7.65/17) 


45+0-09 


< 0.04(0.00) 


0.47 (7.93/17 






LMC 


-0 31+0-80 


< 0.35(0.11) 


0.38 (6.52/17) 


45+0-09 


< 0.02(0.00) 


0.47 (7.93/17 






SMC 


04+0-45 


< 0.15(0.04) 


0.40 (6.88/17) 


45+0.09 


< 0.02(0.00) 


0.47 (7.93/17 



GRB 061121 


50, 1300 


MW 


0.72i0.2J 


< 0.10(0.00) 


0.96 (15.34/16) 


47+0.00 

'-'•^'-0.12 


< 0.13(0.01) 


1.11 (17.75/11 






LMC 




< 0.69(0.18) 


0.93 (14.82/16) 


•^5+0-12 


06+0-03 

U.UD_g Q5 


0.94 (15.06/11 






SMC 


-i-i5ll;Il 


19+0-14 


0.72 (11.55/16) 


U-JO_g gg 


04+0-02 


0.87 (13.96/11 



GRB 061222A 


126 


MW 
LMC 
SMC 


95+0-07 
95+0-07 
95+0-07 


90+0-30 

90+0-4" 

1 00+"-'"' 

1-^^-0. 30 




GRB 070306 


1400, 5800 


MW 
LMC 
SMC 


39+0-13 
39+0-13 

U-Oy_Q gg 
„ OQ + 0.13 

U-59_o gg 


1 Q5+0.52 

J--oo_0.46 
no+0.52 

^■U^_g 51 


OO (1.40/0) 
OO (1.40/0) 
OO (1.40/0) 


GRB 070328 




MW 









LMC 
SMC 



GRB 070521 


90 


MW 
LMC 
SMC 








0.40 
0.40 
0.40 


> 0.8 

> 0.7 

> 0.8 




GRB 071020 


72, 220 


MW 0.58l?;28 

LMC 

SMC 


< 0.26(0.00) 


1.40 (20.97/15) 

< 0.36(0.00) 

< 0.40(0.00) 


75+0-20 
1.40 (20.97/15) 
1.40 (20.97/15) 


< 0.04(0.00) 

„ 75+0.20 
f, 75+0.20 

'0_g.gg 


1.44 (21.54/15) 

< 0.05(0.00) 

< 0.05(0.00) 


1.44 (21.54/lJ 
1.44 (21.54/11 


GRB071112C 


5.8,22 


MW 
LMC 
SMC 




< 0.17(0.00) 

< 0.13(0.00) 

< 0.10(0.00) 


1.03 (33.86/33) 
1.03 (.33.86/33) 
1.03 (33.86/33) 


49+0-01 

49+0-01 

49+0-01 
"•^"-0.37 


< 0.17(0.00) 

< 0.13(0.00) 

< 0.10(0.00) 


1.03 (33.86/3: 
1.03 (33.86/3: 
1.03 (33.86/3: 



GRB 071117 



MW 

LMC 
SMC 



GRB 080319B 


72, 860 


MW 


0.39iro^ 


< 0.01(0.00) 


1.12 (110.87/99) 


n oo+o.oo 


< 0.01(0.00) 


1.12 (110.98/9 






LMC 


0.38iro^ 


< 0.02(0.00) 


1.12 (110.87/99) 


n oo+o.oo 


< 0.02(0.00) 


1.12 (110.98/9 






SMC 


0.221^15 


03+0-02 

'-'■^•^-0.02 


1.09 (108.32/99) 


26+0-11 


< 0.02(0.02) 


1.10 (108.46/9 



GRB 080319C 


18, 36 


MW 


1 61+0-02 


06+0-05 


0.89 (1.79/2) 


0.75 


0.13 


2.52 (5.04/2; 






LMC 


1 07+0.86 


11+0-09 


0.89 (1.79/2) 


0.75l0g-00 


-0.02 


1.41 (2.82/2; 






SMC 


0A8t\il 


26+0-23 


0.89 (1.79/2) 


0.4810-27 


26+0-00 

"■^°-0.07 


0.89 (1.79/2; 



GRB 080413B 


19,100 


MW 


„ 2S+0.04 


02+0 01 
"•"^-0.01 


0.92 (140.08/153) 


40+0-01 
'-'•^'-'-0.00 


< 0.01(0.00) 


1.09 (166.49/lJ 






LMC 


19+0-07 


03+0-01 


0.91 (139.90/153) 


40+0-01 


< 0.01(0.00) 


1.09 (166.49/lJ 






SMC 


n ,0+0.08 
^■'-■^-0.08 


04+0-01 
L'-"''_0.01 


0.92 (140.28/153) 


40+0-01 


< 0.01(0.00) 


1.09 (166.49/11 



GRB 080430 



1.4, 180 MW 
LMC 
SMC 



0.69 
0.50 
1.12 



+0.60 
-0.66 
+0.95 
-1.27 
+0.37 
1.00 



12+0-12 

'-'-1^-0.11 

< 0.28(0.10) 

< 0.15(0.02) 



1.30 (5.19/4) 
1.65 (6.62/4) 
1.86 (7.45/4) 



50+0-12 
^-0^^-0.06 

49+0-13 
"•^^-0.00 

0.62 



14+0-04 

"■1^-0.06 
10+0-02 

0.06 



1.68 (6.72/4; 
1.65 (6.62/4; 
1.98 (7.94/4; 



GRB 080602 MW 

LMC 
SMC 



GRB 080603B 


160, 720 


MW 
LMC 
SMC 


52+0-25 

"■3^-0.30 
52+0-25 

n 00+0.55 
"■^^-0.74 


< 0.06(0.00) 

< 0.35(0.00) 

< 0.16(0.05) 


0.41 (1.22/3) 
0.41 (1.22/3) 
0.26 (0.77/3) 


54+0-09 

54+0-09 

19+0-44 
"•^^-0.03 


< 0.06(0.00) 

< 0.20(0.00) 

< 0.06(0.05) 


0.41 (1.23/3; 
0.41 (1.23/3; 
0.26 (0.77/3; 


GRB 080605 


72,220 


MW 
LMC 
SMC 


1.57 

0.98 

65+0-24 


0.00 

0.12 

16+004 


3.42 (30.81/9) 

2.52 (22.69/9) 
0.28 (2.49/9) 


0.97 
0.97 

74+0.20 

"•'^-0.04 


0.12 
0.12 
14+0-01 


2.86 (88.34/1: 
2.52 (22.69/9 
0.30 (2.73/9; 
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Table AlO - continued 



GRB 080607 


13, 22 


MW 
LMC 
SMC 


-1 7O+0.27 
-^■'''-0.29 

1.28 
2.94 


-^0+0-09 
u.au_Q (j7 

0.31 

0.00 


0.87 (7.81/9) 
4.50 (40.49/9) 
5.82 (52.37/9) 


1 ig+o.oo 

1.19 
1.19 


43+0.03 
"-^'-'-0.02 

0.33 
0.25 


1.46 (13.14/9) 
4.50 (40.53/9) 
11.96 (107.67/9) 


GRB 080613B 




MW 
LMC 
SMC 














^jrlXJj UoU 1 Z 1 


OoU, OoUUU 


IVl vv 

LMC 
SMC 


_4 13+5.13 


< 0.60(0.21) 

< 1.05(0.42) 


n A'\ ("^ on /Q^ 
0.42 (3.74/9) 
0.39 (3.48/9) 


n sfi+o-io 
0-80-0. 00 
n sfi+o.io 
U-»o_o.oo 


< 0.26(0.07) 

< 0.12(0.01) 


n APi { A ^fi /Q^ 

u.^D y^.io/ y ) 
0.52 (4.65/9) 
0.57 (5.10/9) 




0.C5, iZ 


iVl VV 

LMC 
SMC 


_-, J. +3.39 
'^-10.08 


^ U-Ooi^vj.uy ) 

< 1.51(0.25) 

< 1.76(0.59) 


0.49 (1.97/4) 
0.49 (1.97/4) 


n 3R+0.23 

n 07+0.22 
"••^'-0.02 


f, nq+o-os 
0-uy_o.o5 

o.oslH^ 

06+0-03 

"■O"-0.04 


0.51 (2.03/4) 
0.52 (2.07/4) 


R R nsfiQi ^^ A 


0, zyu 


iVl VV 

LMC 

SMC 


-, 47+0.09 

1 47+0.09 

1 47+0.09 
-^•^'-0.43 


<^ U-ZZl^U.UU ) 

< 0.26(0.00) 

< 0.26(0.00) 


1.06 (7.44/7) 
1.06 (7.44/7) 


-I ,p,+0.00 
^•^0_„ 26 

1 20+0-00 
1 20+0-00 


n „7+0.08 
0-U'-o.02 

07+0-00 
^'-^'-0.02 

07+0-00 
^-o' -0.02 


±.ou t^y.u 1 / i ) 
1.24 (8.65/7) 
1.24 (8.67/7) 


kjrtJj UoiUU i 


in 1 9(1 
lU, IZU 


iVl VV 

LMC 
SMC 


A qf; + 0.28 
^■•^0-0.33 
31+0.37 

43+0-33 


n -13+0.09 
^•-^■'-0.07 

11+0-09 


n QR /''^'^ '?S/'^7^ 
u.yo yOO.OOlOt ) 

0.99 (36.68/37) 
1.04 (38.49/37) 


n «fi+o-06 

U.SD_o Qo 

n sfi+0-07 
U.86_o Qo 

n sfi+O-07 


< 0.04(0.01) 

< 0.04(0.00) 


^ ^'i ( A^ q7 I'\7 
±.10 ^^^i.y (lot 

1.13 (41.97/37 

1.14 (42.15/37) 


GRB 081121 


43, 430 


MW 
LMC 
SMC 


_n 1 3+0.25 
'^■^'^-0.28 
-0 38+0-50 
-0 31+0-43 


n in+o-03 

'J-^'J-0.03 

n 15+0-05 

11+0-03 


0.65 (3.93/6) 
0.29 (1.74/6) 
0.53 (3.18/6) 


n 37+0.13 
f, 37+0.11 

o.37l°;Jg 


0.08l0;0i 

07+0-01 
-0.01 

06+0-01 
u.uo_o 01 


0.95 (5.71/6) 
1.43 (8.55/6) 
1.48 (8.87/6) 


GRB 081203A 


140, 7200 


MW 
LMC 
SMC 




0.1710-0^ 
0.29+0-i« 
49+0-24 


0.75 (5.27/7) 
0.74 (5.15/7) 
0.74 (5.15/7) 


0-66ir,;i^ 

57+0-02 
"•°'-0.17 

n cv+0.12 

0.671007 


n -17+0.03 
n -15+0-02 

^'■l''-0.03 
12+0-02 
'^-1^-0.02 


0.75 (5.27/7) 
0.92 (6.42/7) 
1.16 (8.14/7) 


GRB 081221 


190 


MW 
LMC 
SMC 








0.89 
0.89 
0.89 


> 0.5 

> 0.5 

> 0.5 




GRB 081222 


1.2, 43 


MW 

LMC 
SMC 


-0.45 

-0 31+0.70 
-0 01+0-52 


0.02 

15+0-12 
'J-l''-0.12 

< 0.19(0.07) 


1.45 (27.54/19) 
1.30 (24.68/19) 
1.33 (25.19/19) 


0.48 

47+0-12 
^•^'-0.00 

47+0-12 
-0.00 


< 0.03(0.02) 

< 0.03(0.01) 

< 0.02(0.01) 


1.45 (27.54/19) 

1.43 (27.19/19) 
1.42 (27.03/19) 


GRB 090102 


4.3, 29 


MW 

LMC 
SMC 


1 X4 

U.90 Q gr^ 

-0 16+^-00 

"■'-"-1.00 


< 0.20(0.00) 
21+0-21 
11+010 


1.06 (6.33/6) 
0.67 (4.00/6) 
0.54 (3.23/6) 


o-36i;j:;f, 

22+0-1^* 
'^•^^-o.oo 

27+0-09 


n 14+0.06 

'J- 1^-0.03 

14+0-02 

11+0-03 
"■^^-0.03 


1.71 (10.27/6) 
0.71 (4.27/6) 
0.54 (3.23/6) 


GRB 090201 


434 


MW 
LMC 
SMC 








0.49 
0.49 
0.49 


> 0.4 

> 0.4 

> 0.4 




GKB 090424 


120, 4300 


MW 
LMC 
SMC 


-, oK+0.25 
l.ciO_(, 27 

-I qr,+().29 
-^■•^^-0.21 
1 41+0.20 


n nn+0.09 
0-09-0.08 
09+0-06 

06+0-04 

U.UD_(, 04 


0.85 (27.24/32) 
0.75 (23.97/32) 
0.74 (23.80/32) 


n KK+O.OO 
U.00_o 06 

0.55 
0.55 


n oc:+0-02 
u.oo_o 01 

0.31 

0.31 


1.31 (41.81/32) 
1.50 (47.97/32) 
2.15 (68.70/32) 


GRB 090709A 


1.4, 14 


MW 
LMC 
SMC 


-0.03r7.(i4 
"■^^-7.48 

1.52 


< 5.88 

< 5.02 
0.45 


1.84 (9.18/5) 
1.84 (9.19/5) 
1.87 (9.36/5) 


n or; + 0.09 
'J -^^-0.05 

44+0-00 
^•^^-0.14 

0.44 


n O-1+0.16 
U.»3_o.l6 
75+0-18 

•J-'^-o.ii 
0.77 


1.84 (9.19/5) 
1.84 (9.20/5) 
1.90 (9.49/5) 


GRB 090715B 


8.6, 430 


MW 
LMC 
SMC 


1 40+0-12 
^•^^-0.12 
11+0-38 
"•^^-0.44 
-1 09+0-85 


r, 1O+0.07 

17+0-07 

r, ,0+0.08 


0.67 (31.94/48) 
0.66 (31.84/48) 
0.66 (31.78/48) 


0.63 

45+0-00 
45+0-00 


0.00 
11+0-01 
06+0-01 


2.24 (107.63/48) 
0.69 (33.20/48) 
0.79 (38.07/48) 


GRB 090812 


1.4, 290 


MW 
LMC 
SMC 


1 24+0-2S 
0.581";- 


< 0.16(0.03) 

< 0.33(0.12) 

< 0.24(0.08) 


1.52 (12.15/8) 
1.31 (10.51/8) 
1.38 (11.08/8) 


0.52 

51+0-01 
'J-Ol-0.12 
52+0-00 


0.16 

0-121;:;;]^ 
0.08l0;03 


2.65 (21.20/8) 
1.31 (10.52/8) 
1.39 (11.09/8) 


GRB 090926B 




MW 
LMC 
SMC 
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Table AlO - continued 







MW 
LMC 
SMC 


57+0-"! 

,.-+0.01 


< 0.01(0.00) 

< 0.02(0.00) 


1 1 5 ("QS 54/81 
1.15 (93.51/81) 
1.15 (93.51/81) 


57+0-01 
'^•>''-0.02 
57+0-01 

^7+0.01 


<- 01 ("0 no~i 

< 0.01(0.00) 

< 0.02(0.00) 


1 16 ("93 64/81 
1.16 (93.64/81) 
1.16 (93.64/81) 




12, 720 


MW 
LMC 
SMC 


"^■"^-0.16 
n oQ+0.65 

2.38 


10+°°2 
^J-l^"_0.03 

n ■J5+0.10 
0.00 


1 40 (^23 7Q /1 7") 
1.39 (23.61/17) 
2.70 (45.98/17) 


0.66 

n 55+0-11 

"J-OO-O.OO 
0.66 


0.29 

oo+O-Ol 

0.18 


1 3 QQ (^237 76/1 7^1 
1.42 (24.12/17) 
4.45 (75.62/17) 


GRB 091127 


86, 240 


MW 

LMC 
SMC 


r. oQ + 0.02 

^•^•^-0.02 

n r,o + 0.03 

"•^•^-0.03 
n 9S+0.03 
U-^»-0.02 


n n/+0.()l 

04+° °^ 

U-U^.Q.Ol 

02+° °^ 
'^■^^-0.01 


0.93 (367.11/396) 
0.95 (374.56/396) 
0.98 (387.70/396) 


n on + 0.02 

22+°-04 

U-ZZ_Q Q2 

29+0-00 
^•^^-0.04 


n n/i+O.Ol 

04+0-01 
02+0-01 


0.93 (367.23/396) 
0.95 (375.96/396) 
0.98 (387.95/396) 


GRB 091208B 


580, 1200 


MW 
LMC 

SMC 


-^•^^-0.64 
-^•US_0.88 


< 0.19(0.07) 

< 0.29(0.10) 

< 0.43(0.14) 


1.08 (6.51/6) 

1.08 (6.48/6) 

1.09 (6.54/6) 


1 07^ 

1 Q7+o.()() 
94+0-13 


0-091^3 

"--^^ -0.03 

14+0-06 

"■-^^-0.07 


1.11 (6.66/6) 

1.08 (6.48/6) 

1.09 (6.54/6) 


GRB 100615A 


24 


MW 
LMC 
SMC 








0.76 
0.76 
0.76 


> 2.8 

> 2.8 

> 3.2 




GRB 100621 A 


100, 140 


MW 
LMC 
SMC 


1.07+071 

1 19+0-64 

1.80 


0.87 


1.77 (7.08/4) 
1.84 (7.37/4) 
2.52 (10.09/4) 


-1 nn + O.Ol 

1.02+,, j2 
i.Uci_Q 25 

1.03 


1 ,„+0.12 
l-l'-().04 

1 10+0-11 

l-l'J_o 03 

1.16 


1.77 (7.09/4) 
1.88 (7.52/4 
3.49 (13.94/4) 


^jrrxo lUU 1 ZOO 


9 9 1 An 


IVi vv 

LMC 
SMC 


n /ic+0.50 

n OQ+0.66 
"•^»_0.81 

70+°-3i 


< 0.35(0.10) 

< 0.17(0.02) 


yz. 1 Kj/ O) 
0.45 (2.68/6) 
0.60 (3.58/6) 


o 47+0.28 
^•^'-0.07 

"•^u_0.00 
70+0-05 


< 0.14(0.08) 

< 0.11(0.02) 


yZ. 1 Kj/ D) 

0.46 (2.74/6) 
0.60 (3.58/6) 


GRB 110205 A 


360, 580 


MW 
LMC 
SMC 


1.23 
55+0-10 

0.78 


0.00 

"•^^-0.01 

0.04 


1.72 (139.47/81) 
1.18 (95.39/81) 
1.28 (104.01/81) 


0.72 

55+0-14 
0.72 


0.12 
11+0-01 

"•^J^_0.02 

0.05 


9.00 (729.07/81) 
1.18 (95.40/81) 
1.29 (104.57/81) 


GRB 110503A 


4.3, 720 


MW 
LMC 
SMC 


1 24+0-4" 

1 00+0.36 
l-^»_0.74 


< 0.23(0.02) 

< 0.24(0.00) 

< 0.14(0.00) 


1.27 (2.55/2) 
1.31 (2.63/2) 
1.31 (2.63/2) 


0.49 

49+0-00 
0.49 


0.17 
10+0-07 

0.06 


4.64 (9.27/2) 
2.21 (4.42/2) 
2.45 (4.90/2) 
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